On Collisionless Electron-Ion Temperature Equilibration in the 

Fast Solar Wind 

J. Martin Laming 1 



m 
o 
o 

(N 



ABSTRACT 



O 



We explore a mechanism, entirely new to the fast solar wind, of electron heat- 
ing by lower hybrid waves to explain the shift to higher charge states observed in 
Q ■ various elements in the fast wind at 1 A.U. relative to the original coronal hole 

1/-) ■ plasma. This process is a variation on that previously discussed for two temper- 

ature accretion flows by Begelman & Chiueh. Lower hybrid waves are generated 
by gyrating minor ions (mainly a-particles) and become significant once strong 
^ ion cyclotron heating sets in beyond 1.5_R Q . In this way the model avoids con- 

flict with SUMER electron temperature diagnostic measurements between 1 and 
1.5i? . The principal requirement for such a process to work is the existence 
of density gradients in the fast solar wind, with scale length of similar order to 
the proton inertial length. Similar size structures have previously been inferred 
r£3 ■ by other authors from radio scintillation observations and considerations of ion 

cyclotron wave generation by global resonant MHD waves. 

O 

& 
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1. Introduction 

The 1990 October launch of Ulysses and the subsequent polar passes in late 1994, 1995, 
2000 and 2001 have provided a wealth of new data and insights concerning solar polar coronal 
holes and the fast solar wind emanating from them. The south polar pass of 1994 (during 
solar minimum) highlighted a particular problem concerning the charge states observed in 
various elements in the fast wind, which are generally characteristic of higher temperatures 
than observed spectroscopically in coronal holes. Coronal holes show very little emission 
in spectral lines corresponding to plasma temperatures of 10 6 K or greater (e.g. Doschek 
& Feldman 1977; David et al. 1998; Wilhelm et al. 1998), whereas Geiss et al. (1995) find 
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freeze-in temperatures from the O, Si and Fe charge states in the range 1.2 — 1.5 x 10 6 
K. Fe in the fast solar wind is dominated by charge states Fe 10+ and Fe 11+ , whereas lines 
from Fe X, XI, and XII are generally not detectable in spectroscopy of coronal holes; an old 
result going back to Doschek & Feldman (1977). Geiss et al. (1995) and Ko et al. (1997) 
recognized that the electron temperature must increase outwards as the fast wind flows out 
of the coronal hole, and were able to derive a empirical temperature profile, but did not 
specify a mechanism to produce such an effect. 

The heuristic temperature profile of Ko et al. (1997) predicts a maximum in the elec- 
tron temperature of about 1.5 x 10 6 K relatively close to the solar surface, at about 1.5i? 
heliocentric distance. However observations in a temperature diagnostic line ratio in Mg IX 
(Wilhelm et al. 1998) out to l.QR & with the SUMER instrument on SOHO failed to detect 
such a temperature increase. Consequently, in a series of papers (Esser & Edgar 2000, 2001; 
Chen, Esser, & Hu 2003), the effects of a halo electron distribution (i.e. a power law tail 
on the core Maxwellian) that remains consistent with SUMER electron temperature mea- 
surements and differential flows between ions of the same element were investigated. While 
carefully chosen differential flows can explain most of the in situ charge state observations 
without the need for extra electron heating, the required speed differentials are so large as to 
be implausible, leaving electron heating as the most likely explanation. The fast solar wind is 
often observed (at 1 AU) with T e{{ > T e± (Feldman et al. 1975) and with T e \\/T e± > T p \\/T p± , 
and this electron temperature anisotropy can be described in term of a core Maxwellian 
distribution and a halo higher temperature Maxwellian or power law distribution, which 
contains typically 5% of the electrons (Marsch 1991). It has been suggested that this halo 
distribution persists all the way back to the sun (Esser & Edgar 2000), and this might 
be appealing since the electron temperature diagnostic used by Wilhelm et al. (1998) (Mg 
IX) belongs to an isoelectronic sequence that is in general insensitive to such a population 
(Keenan 1984), whereas the ionization balance is not. However the O VI diagnostic used 
by David et al. (1998) should be sensitive to suprathermal electrons but shows no evidence 
of their existence out to 1.3 R Q . In any case, it seems that such an explanation must fail, 
because the halo distribution will transfer energy to the core Maxwellian by Coulomb col- 
lisions on a timescale ~ (10 7 — 10 s ) x T e 3 / 2 /n e s, (T e6 is the electron temperature in 10 6 K, 
and n e is the electron density in cm -3 ) and even with the fastest inferred outflow speeds, 
the solar wind would not get very far before a discernable increase in the core Maxwellian 
electron temperature should become apparent. 
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2. Wave Acceleration and Heating in the Fast Solar Wind 

The acceleration of the fast solar wind is widely believed to occur through ion cyclotron 
waves, though the manner in which they are excited remains controversial. The original 
concept, that low frequency waves alone are sufficient to accelerate the fast solar wind to its 
eventual speed of 700-800 km s _1 , appears to be adequate for acceleration of protons but not 
for the minor ions (Ofman & Davila 2001), which in this case would be accelerated mainly by 
the Coulomb drag of the protons. Hence minor ions should flow out of the coronal hole with 
lower speeds than the protons. Renewed indications that ion cyclotron waves are required 
come from the UVCS Doppler dimming measurements of the 5+ outflow velocity, which 
are large compared to hydrogen, ~ 300 — 400 km s _1 at 2.5-3.0 i? heliocentric distance. By 
contrast a purely low frequency wave driven outflow would produce 5+ outflow velocities 
of less than 100 km s -1 at this distance. Similar indications come from the observation 
in the fast solar wind of faster flow speeds for a particles (and in fact all minor ions) 
than for protons with Ulysses and Helios (Reisenfeld et al. 2001; von Steiger et al. 2000; 
Neugebauer et al. 1996). The Ulysses and Helios observations only go in to the sun as far 
as ~ \ARq heliocentric distance, but the inference that a particles flow out from the sun 
at least as fast as protons in the fast solar wind may be extended to lower altitudes by 
spectroscopic observations with SUMER of the He abundance (Laming & Feldman 2003). 
Here He abundances at or below the usual fast solar wind value of 5% relative to H are found, 
implying flow speeds at these lower altitudes at or above the H flow speed to keep the flux 
of a particles constant. This is also consistent with recent modelling by Li (2003), where by 
heating primarily the a particles, a steep transition region, hot corona and fast wind can be 
created by high frequency Alfven waves propagating up from the coronal base. On the other 
hand Cranmer, Field, & Kohl (1999) and Cranmer (2000) argue that the high frequency 
(or ion cyclotron) waves that accelerate the fast wind must be generated throughout the 
extended corona, due to their rapid dissipation once generated, and the observation that 
ions continue to be heated in the range 1.5 — 5R Q heliocentric distance. 

We propose an entirely new (to the fast wind; though see Schwartz, Feldman, & Gary 
1981, for related considerations) means of heating/accelerating electrons. At a heliocentric 
distance between 1.5 and 2 R & , the perpendicular velocities of the minor ions start to exceed 
those of the protons (e.g. Cranmer et al. 1999). There are also indications from observations 
that a particles are similarly heated. All these ions will have significantly larger gyroradii 
than the protons. In the presence of a density gradient, these ions may then excite lower 
hybrid waves in the colder protons. This requires that the electron gyroradius be less than 
the wavelength/27r. The waves will damp by heating electrons in a direction parallel to 
the magnetic field. Unless other restrictions set in, the waves ultimately will saturate once 
the electrons are heated such that their gyroradii are no longer sufficiently small. In the 
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region of interest the plasma conditions are approximately magnetic field B — 1 G and 
density n e = 10 6 cm~ 3 . Then the proton and electron gyrofrequencies are Q p = 10 4 and 
f2 e = f .8 x 10 7 , and the lower hybrid frequency is VL LH = ^Vl p VL e = 4 x 10 5 (all in rad s _1 ). 
The electron and proton plasma frequencies are ui pe = 6 x 10 7 and oo pp = 1.5 x 10 6 . It will turn 
out that waves are excited at low k by ions moving with velocities greater than the proton 
thermal speed, so the limit uj » \/2kvi of the plasma dispersion function is appropriate, 
and that a particles will be species most important for the wave generation. This is due 
to their large gyroradius compared to protons and their high abundance compared to other 
minor ions. 

Lower hybrid waves are electrostatic ion oscillations, which can occur in magnetic fields 
strong enough that the electron gyroradius is smaller than the lower hybrid wavelength/27r. 
Due to this necessary magnetization of the electrons, the waves propagate preferentially 
across magnetic field lines. The parallel component of the wavevector k\\/k < u pi /u pe . Since 
co/k± ~ {m e /rrii) l l 2 u /k\\ the wave can simultaneously be in resonance with ions moving 
across the magnetic field and electrons moving along magnetic field lines, which facilitates 
collisionless energy exchange between ions and electrons on timescales much faster than 
that associated with Coulomb collisions. They have also been discussed in connection with 
cometary X-ray emission, electron acceleration in solar flares, supernova remnant shock waves 
and Advection Dominated Accretion Flows (ADAFS) (Vaisberg et al. 1983; Kranosel'skikh et 
al. 1985; Bingham et al. 1997; Shapiro et al. 1999; Bingham et al. 2000; Begelman & Chiueh 
1988; Luo, Wei, & Feng 2003; McClements et al. 1997), and observed in situ together with 
accelerated electrons at Halley's comet (Gringauz et al. 1986; Klimov et al. 1986). A cold 
plasma theory for lower-hybrid waves is given in the Appendix of Laming (2001a). Here we 
summarize the theory with finite electron and ion temperatures. 

The dispersion relation is (see e.g Laming 2001b) 

co 2 pe (lk 2 /k 2 )/(l + ^Jc 2 k 2 ) 
UJ — 2 — 7 7 \\ d (1) 

1 + &( 1 -*(;*r)) + ifc< 1 - i ) 

where I = J + °° J 2 (^^) exp (^^p) v ± dv ± , v 2 = k B T e /m e and v 2 = k^m^ and 
(z) = —z/^/lr exp (— t 2 ) I (t — z) dt is the usual plasma dispersion function (T e j and 
m e j are electron and ion temperatures and masses respectively and k-Q is Boltzmann's con- 
stant). Specializing to uo » \plkvi and Vl e >> k±v e ±, so that 7 ~ 1 — k\v 2 e J^ e and 
~ 1 + i-sj it/2 (oj/kvi) exp (— oo 2 /2k 2 v 2 ) and taking u pe » f2 e and c — > oo, 

k\\ \ Q 2 e i f^ UJ^l 2 e UJ 2 pi . 2 2 2 x 



(2) 
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The instability is illustrated schematically in Figure I. The density is increasing towards 
the right, and the density gradient is perpendicular to the magnetic field. A local anisotropy 
in the ion distribution develops in the direction formed by the cross product of density 
gradient and magnetic field, which in the representation of Figure 1 is into or out of the page. 
This is similar to a scenario envisaged by Begelman & Chiueh (1988), who were interested 
in electron-ion equilibration in two temperature accretion flows. We consider the effect of a 
density gradient in a uniform magnetic field. In such a situation the ion distribution function 
is given by 

— * 

where the magnetic field B = Bz, and the density gradient {dn/dx)x is related to L by 
n(dx/dn) = L. We consider the growth of lower hybrid waves in the protons excited by 
minor ions with large gyroradii. The growth rate for each minor ion species is given by 
(Laming 2001a) 

7* = & (l + f cos 2 9) " J k- §P (u - k- v) 

exp (-»£-) , (4) 



pi 
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where 6 is the angle between B and k, v thl _ = \Jk^Ti^jmi and Afq 2 /M is the product of 
the element abundance, ionization fraction and charge squared of the particular ions that 
excite the wave, divided by their mass in atomic mass units. Protons are expected to have 
insufficiently large gyroradii to excite lower-hybrid waves, and so the factor Afq 2 /M favors 
a-particles over other minor ions in the wave generation. The term —1 in the square brackets 
represents the ion Landau damping given by the imaginary part of the expression for u in 
equation 2. Equation 4 is similar to equation (3.20) in Begelman & Chiueh (1988), with the 
identification of their drift velocity Vd% = —v 2 h± /QiL, and in their case uj/k « v t h±- Their 
treatment yields waves at much higher k than ours does, since they consider a single ion 
distribution drifting with velocity v di producing waves at k ~ oo/vdi, whereas we treat the 
effect at a single point of a continuum of different Maxwellians with differing densities (but 
otherwise the same) producing waves at k ~ (ou/vth±) i r g/L). This difference is appropriate 
in view of the fact that they are interested in plasmas with ion pressure >> magnetic pressure 
>> electron pressure where the plasma itself amplifies the magnetic field, while here we are 
concerned with magnetic pressure >> ion and electron pressures where the magnetic field 
is imposed externally. Begelman & Chiueh (1988) also considered the effects of magnetic 
curvature. We assume that the necessary gradients in magnetic field are much less likely 
to exist in the low (5 plasma of the solar coronal hole. We find the wavevector k max where 
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the growth rate is maximized in the direction perpendicular to B (i.e. 9 = 7r/2), and the 
maximum growth rate itself, which is plotted in units of w x Afq 2 /M in Figure 2. 

If the waves are driven to marginal stability k max ~ Q e /v e ± ~ ur g /Lv i± , where Vi± is 
the perpendicular ion velocity, equation 1 can be written 

V e± = 7 x " T- ( 5 ) 

^V(^A 2 )+<(i + </^ 2 ) 

At saturation where /c_i_t> e _i_/f2 e ~ 1, J ~ 1/2, and putting fc||/A; = cos#, 

2 _ jj^ (1 + tjg/gg) (1 + ^Agg) L 2 
Ue± cos 2 0) /2 + < (1 + u&t&A^) r| ■ U 

This is a quadratic equation for v 2 ± which has lowest order solution 

„ wleVf*. ( UJ 2 „ \ _1 L 2 

The kinetic growth rate varies as 7 oc f 1 + cos 2 0^ , (Laming 2001a) and averaging v 2 ± 

V w pp / ■ 

over this growth rate gives (v 2 ± ) = v 2 ± (l — 1/V2) tu 2 e L 2 /tu 2 p r 2 . Thus in principle electron 
and ion temperatures can equilibrate (assuming that (l - 1/V2) L 2 /r 2 > 1), in the sense 
that the waves do not saturate before this is achieved. 



3. Solar Wind Ionization Balance Models 

We use an adaptation of the BLASPHEMER (BLASt Propagation in Highly EMit- 
ting EnviRonment) 2 code (Laming 2001b; Laming & Grun 2002, 2003; Laming & Hwang 
2003), which follows the time dependent ionization balance and temperatures of a Lagrangian 
plasma parcel as it expands in the solar wind. The density rii q of ions of element % with charge 
q is given by 

^ (Cjon,?— q—1 Cion^^i?) ~\~Tl e (CVr,<ji+l H~ C<2r,g+l) 7lj <jr+l Tl e (Crr,q Cdr,q) n>iq (8) 

where Ci oritq , C„ S) Cd r , q are the rates for electron impact ionization, radiative recombination 
and dielectronic recombination respectively, out of the charge state q. These rates are the 



2 The name gives away its origin in modelling laboratory and astrophysical shock waves. 
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same as those used in the recent ionization balance calculations of Mazzotta et al. (1998), 
using subroutines kindly supplied by Dr P. Mazzotta (private communication 2000). The 
electron density n e is determined from the condition that the plasma be electrically neutral. 
The ion and electron temperatures, T iq and T e are coupled by Coulomb collisions by 

and 

dT e _ 0.13n e {T iq - T e ) q 2 n iq / (g + 1) / In A \ T e / dn e \ 2 dQ ^A lig U w 

dt T ? /2 X M ^ (V,,, //,,,) V 37 n e \dt) ion 3n e k B dt ^ 3 n e k B 

(10) 

Here M iq is the atomic mass of the ions of element % and charge q in the plasma, and In A ~ 28 
is the Coulomb logarithm. The term in dQ/dT represents plasma energy losses due to ion- 
ization and radiation. Radiation losses can be taken from Summers & McWhirter (1979), 
though are generally negligible in applications to the solar wind. The term — [T e j n e ) (dn e / dt) ic 
gives the reduction in electron temperature when the electron density increases due to ion- 
ization. Recombinations, which reduce the electron density, do not result in an increase in 
the electron temperature in low density plasmas, since the energy of the recombined electron 
is radiated away (in either radiative or dielectronic recombination), rather than being shared 
with the other plasma electrons as would be the case for three-body recombination in dense 
plasmas. 

The last terms in equations 9 and 10 represent collisionless ion-electron energy transfer, 
which we estimate assuming that the population of lower hybrid waves is approximately 
steady state. Then the wave growth rate due to the gyrating ions is equal to the wave 
Landau damping rate due to the electrons, and the energy transfer rate is given by 2^iU w 
in ergs cm~ 3 s _1 where U w = Eq/8tt is the wave energy density in terms of the wave peak 
electric field Eq. A threshold electric field for lower hybrid waves between the linear and 
stochastic regimes has been determined by Karney (1978) by numerical integration of the 
Hamiltonian equations of motion, and also analytically using an ion trapping argument, to 
be 



4 y uj ) kc 

where B Q is the background magnetic field and fl p is the proton gyrofrequency. This field 
marks the onset of proton heating as well as electron heating by the waves. We take E = 
E t h r (probably a conservative assumption). Writing the electron heating rate in terms of 
both collisional and collisionless processes, with >> T e and dropping the q subscript in 
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specializing to one ion species alone, 




+ 



n e 647r \ u> J rriiC 




(12) 



where k max is the wavevector where the growth rate is maximized, and Tj = 2Tjj_/3 = 
rriivf ± /3 = 2mivf y /Z since v? ± = v\ x + vf y = 2vf y » v^. In our calculations we modify In A 
from its usual value to account for the collisionless energy transfer between a-particles and 
electrons, its new value being given by 



the quantity ^iMi/toAfq 2 being plotted in Figure 2. The same value is applied to electron 
heating by other minor ions, and the original value of In A is kept for proton-electron energy 
transfer. In the region of the solar wind around 1.5R@ heliocentric distance, where T e ~ 10 6 
K, B ~ 1 G, n e ~ 10 6 cm -3 and 7' = , ~fiM i /ujAfq 2 (u/kvi y ) 2 ~ 0.1 — 1 for a-particles of 
around 5% abundance relative to protons, the electron heating rate is increased typically by 
factors of order 100 - 1000 over that due solely to Coulomb collisions with ions heavier than 



We use the analytic model of Banaszkiewicz, Axford & McKenzie (1998) for the magnetic 
field strength. We take the perpendicular velocity for all ions heavier than H from the 
empirical model for O VI (model B2; equation 28) in Cranmer et al. (1999) in the range 
1.5i? — 3.5i? . Below this we interpolate between values of 20 km s _1 on the solar surface 
and 84 km s _1 at 1.5i? . The flow speed is linearly interpolated between its initial value at 
the solar surface and 130 km s _1 at 1.7R Q , a value determined by Giordano et al. (2000). 
Above 1.7 Rq the flow speed is allowed to evolve according the the action of the adiabatic 
invariant in the diverging magnetic field with the perpendicular velocity specified as above, 
and above 3.5i? the perpendicular velocity is also allowed to evolve in this manner, rather 
than being determined by the fits in Cranmer et al. (1999). The resulting flow speed is a 
very good match to the empirical models Bl and B2 in Cranmer et al. (1999), as shown in 
Figure 3 for various values of the initial flow speed, determined to be in the range 3-60 
km s' 1 by various authors (Patsourakos & Vial 2000; Hassler et al. 1999; Wilhelm et al. 
2000; Gabriel, Bely-Dubau, & Lemaire 2003). Figure 4 shows the electron density in the 
simulations for initial flow speeds between 5 and 60 km s -1 , compared with measurements 
using a diagnostic line ratio in Si VIII by Doschek et al. (1997) and from UVCS polarization 
brightness observations given by Cranmer et al. (1999). An initial flow speed of around 10-20 
km s _1 appears to be the best match to the various observations. 

Equations 8-10 are integrated following the solar wind out from an initial position at 
l.O5i? out to around 7R G by which time charge states are frozen in. Each model follows H, 




(13) 



H. 
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He, and one minor ion, C, O, Mg, Si, or Fe in the present work, with fractional abundances 
0.83, 0.16 and 0.01 respectively by mass. The initial electron and ion temperatures are 
9 x 10 5 K, which also establish the initial ionization balance. The density at this point is 
taken to be 10 8 cm~ 3 . After each time step in the ionization balance, the densities and 
electron temperature are modified according to the adiabatic expansion of the solar wind 
governed by the magnetic field geometry and the wind velocity as specified above. 

Figure 5 shows the electron temperature profile resulting from simulations with an 
initial flow speed of 10 km s _1 and various values of the lower hybrid growth rate, which 
depends on the density gradient assumed. For reference, the temperature measurements 
of Wilhelm et al. (1998) using the Mg IX 750 A /706A diagnostic line ratio are also given. 
The error bar on the points at 1.3 and 1.6i? are estimated here from the scatter in the 
points on Figure 8 of Wilhelm et al. (1998). These authors use the atomic physics data of 
Keenan et al. (1984) to derive a temperature ratio, which is consistent with coronal hole 
temperatures determined by other authors (e.g. Doschek & Feldman 1977; Doschek et al. 
2001). More recent calculations in R-matrix and distorted wave approximations summarized 
by Landi et al. (2001) give temperatures significantly higher or lower respectively, and are 
likely due to inaccuracies in these atomic data. We base our coronal hole temperatures on 
works such as Doschek & Feldman (1977) and Doschek et al. (2001) where the ionization 
balance is the principal temperature diagnostic, and on the O VI observations of David et 
al. (1998), whose electron temperature diagnostic depends on much less controversial atomic 
physics. Even if the absolute temperatures measured by Wilhelm et al. (1998) cannot be 
interpreted with confidence, as they state in their paper, at least the maximum variation of 
the electron temperature with distance from the solar surface can be constrained by their 
observations. Figure 6 shows the electron temperature profiles resulting from using the 
maximum collisionless energy transfer in Figure 5 and differing initial flow speeds, with the 
highest temperatures resulting from the slowest initial speeds since more time is available 
for the plasma electron to be heated and the rate of temperature decrease due to adiabatic 
expansion is lower. 

Figures 7, 8, and 9 show the evolution of the ionization balances of O, Si, and Fe with 
heliocentric distance. In each case the flow starts at a density of 10 8 electrons cm -3 and a 
temperature of 9 x 10 5 K. The initial flow speed is 20 km s -1 and the electron-ion equilibration 
parameter is 7' = ^iMi/toAfq 2 (to 2 /k 2 v 2 y ) = 0.5. Increased ionization commences at around 
1.5i? , in response to the onset of ion cyclotron heating at this location, and charge states 
freeze in between 2 and 2.5R Q at the values found in situ by Ulysses (Geiss et al. 1995; Ko 
et al. 1997). 
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4. Discussion 

4.1. Charge State Distributions 

The resulting ionization balances are given in Tables 1-4 for C, O, Mg, Si and Fe 
respectively. From the C and O ionization balances in Table 1, it is clear than only for 
values of the parameter 7' = {^/iMi/uiAfq 2 ) (uj/kv iy ) 2 ~ 0.5 and initial wind speeds of 10-20 
does sufficient electron heating occur to bring the modelled charge states into agreement 
with those observed. Tables 2-4 verify that these conclusions do not change when other 
elements Mg, Si, and Fe are considered. We will discuss the significance of 7' below. Before 
doing so we remark on the other salient features of the model. The data of Geiss et al. 
(1995) come from a high speed stream observed closer to the ecliptic plane than the polar 
observations of Ko et al. (1997), and so are slightly less highly ionized as would be expected. 
Even so, both sets of observations for the parameters discussed above agree very well with 
the modelled charge state fractions for C, O, Si and Fe, at least for the two or three charge 
states that dominate the distributions. Mg is still observed to be more highly charged than 
the models predict. However we note that the steady state ionization balance for Mg changes 
dramatically in the temperature range 9 x 10 5 K to 1.1 x 10 6 K (Mazzotta et al. 1998), and 
so we anticipate that small changes in either the model or the atomic data for Mg would 
also bring this element into agreement. We find no need to invoke different outward flow 
speeds for the different elements. C, O, Mg, Si, and Fe all flow at the flow speeds determined 
from O VI Doppler dimming measurements. This is in contrast with the work of Ko et al. 
(1997), where heavier ions needed to flow out at successively lower speeds, as would be the 
case in a thermal conduction driven wind. In addition to a faster outflow, we also model 
higher electron temperatures, as would be expected since less time is available before freeze- 
in, and so higher electron temperatures are required to produce the necessary ionization. 
Even so, the electron temperatures modelled here when extrapolated out to 0.3-1 A.U. are 
still considerably lower than those measured in situ by Helios (Marsch et al. 1989). The 
inclusion of non-zero thermal conduction (see below) might reduce this discrepancy, but it 
would seem that the conclusion of (Marsch et al. 1989) that "electron temperature profiles 
in high-speed streams clearly indicated and even required the existence of heating sources 
other than the one related to the degradation of the electron heat flux" supports our idea 
of electron heating. We also remark that Ko et al. (1997) made some different choices of 
atomic data. In particular their choice of 6+ ionization rate from Lennon et al. (1988) is 
lower than the rate adopted here from Mazzotta et al. (1998) by a factor of about 0.7, and 
elsewhere in the literature (e.g. Moores, Golden, & Sampson 1980; Shull & van Steenberg 
1982) rates for this process can be greater than that of Lennon et al. (1988) by nearly a factor 
of 2. The assessment and validation of atomic data is a huge task, beyond the scope of the 
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current paper (see e.g. Savin & Laming 2002), but clearly central to further quantitative 
development along the lines suggested in this paper. 

Further variations in the atomic rates involved may come from nonthermal electron 
distributions. Throughout this paper we have taken a Maxwellian electron distribution, 
tacitly assuming that the electron-electron collision rate is sufficiently fast to maintain such 
conditions. We argued above that halo electron distributions are unlikely to exist close to 
the Sun because of this collisional equilibration, and numerical estimates suggest that only 
at electron densities < 10 4 — 10 5 cm -3 could a nonthermal electron distribution produced 
by lower hybrid waves survive (by equating the collisional equilibration rate with the wind 
expansion rate). Such densities are only found at or beyond the radial position where ion 
charge states freeze in, and so nonthermal electrons appear unlikely to produce a significant 
change to our charge state results, given the other current observational constraints. However 
lower hybrid waves do appear to be a viable means for producing the electron distributions 
observed in the fast solar wind. 



4.2. Small Scale Structures? 

The determination that 7' = jiMi/uAfq 2 (u/kv iy ) 2 ~ 0.5 — 1 is necessary to produce the 
observed charge states requires the existence of density gradients in the fast wind with scale 
lengths on the order of the a-particle gyroradius, which is about 0.1 km at 1.5 R Q . Radio 
scintillation observations demonstrating the existence of such size scales in the solar wind in 
the ecliptic plane have rather a long history (e.g. Coles & Harmon 1989; Armstrong et al. 
1990; Coles et al. 1991). These are found perpendicular to the magnetic field, with larger 
size scales (typically a factor of 10 within 6 R Q , becoming more isotropic at larger distances) 
inferred along the radial direction. Coles et al. (1995) inferred values of 5n 2 e in polar regions 
at solar minimum to be around 1/10 to 1/15 of that observed in equatorial regions, but due 
to a lack of knowledge of n e in polar regions, were unable to say anything about the variation 
of Sn e /n e . The density measurements reviewed in this paper indicate electron densities in 
polar regions a factor of 1/2 to 1/3 of those in equatorial regions, making 5n e /n e in polar 
regions of similar order to, but still slightly smaller than that in equatorial regions. Absolute 
values of Sn e / n e in coronal hole regions of interest here have been determined observationally 
by Ofman et al. (1997) to be from 0.1 to a few times 0.1. This is smaller than the value ~ 1 
tacitly assumed here, the consequence of which is discussed further below. 

Grail et al. (1997) present more data on the transition from anisotropy inside 5-6 R & on 
scales of order 10 km to isotropy further out, concluding that a real change in the microstruc- 
ture rather than in Alfven wave turbulence takes place, again with reference to the ecliptic 
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plane. Feldman et al. (1996) review these interplanetary scintillation observations together 
with Ulysses observations to constrain the high speed wind structure near its coronal base, 
and argue that the plasma is "sufficiently structured to relax through generation of a drift- 
wave instability that results in electrostatic waves having /c-vectors oriented perpendicular 
to B", which is precisely the motivation for the current work. Grail et al. (1997) go further 
and show that within 6 R Q large scale turbulence is isotropic with a Kolmogorov spectrum 
(structure function oc scale 5//3 ), while smaller scale turbulence shows anisotropy with higher 
structure functions (oc scale) than Kolmogorov turbulence would predict. The scale at which 
this transition takes place which can be interpreted (Woo 1996; Woo & Habbal 1997) as the 
size of the flux tube in which the wind flows, is inferred to be of order 1 km close to the sun. 

In a study of electrostatic ion cyclotron wave generation by global resonant MHD modes, 
Markovskii (2001) suggests that structures may exist with scales down to the proton inertial 
length va/^p = c/u pp , which is of the same order of magnitude as the structures already 
inferred by Woo (1996); Woo & Habbal (1997), and in terms of particle gyroradii about 
an order of magnitude larger than our present inference. However our evaluation of the 
lower hybrid wave growth rate assumed Maxwellian distributions for the minor ions. In the 
presence of a radiation field interacting with the ions, Hasegawa, Mima & Duong-van (1985) 
show that a "kappa" distribution may result; 



/ (v) dv 



n I t /,- • 1 ) r 



2 



1 + 



(14) 



/2W (k-1/2) 

where v t h = ^/k^T/m is the thermal velocity, V is the usual Gamma Function, and k is a 
parameter representing the deviation from a Maxwellian distribution, which is obtained in 
the limit k — > oo. Observations of the a-particle distribution function in the high speed 
solar wind generally give k ~ 3 — 6 (Collier et al. 1996; Chotoo et al. 1998). In Figures 10 
and 11 we plot the dimensionless lower hybrid growth rate •y i Mi/ivAfq 2 and the parameter 
7' = r ) i M i /ujAfq 2 (uj/kviy) 2 against L/r g . From Figure 11 it may be seen the for k — 3 
values of L/r g up to about 7 give acceptable ionization and for lower k even higher L/r g 
would be tolerable. Protons in a kappa distribution may also be able to contribute to the 
wave generation, reducing the requirement on L/r g even more. Given the relatively safe 
assumption that the k value for a-particles in the coronal hole is similar to or less than 
that observed in high speed wind streams at about 1 A.U., we conclude that similar density 
gradients to those proposed by Markovskii (2001) would provide sufficient lower hybrid wave 
generation and electron heating to produce the observed ionization states of minor ions in 
the fast wind. Thus these observations provide further support for the mechanism of ion 
cyclotron wave generation throughout the extended corona proposed by Markovskii (2001), 
that of a global resonant MHD mode driving a cross field current in the resonant layer 
which then excites electrostatic ion cyclotron waves. We have also made rather conservative 
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assumptions concerning the collisionless ion-electron energy transfer. Higher wave electric 
fields than given by equation 11 (up to ~ uB /kc) are certainly possible (Karney 1978; 
Bingham et al. 2003), though as the wave amplitude increases from this point, protons begin 
to be heated as well as electrons (the stochastic regime). Under favorable conditions, the 
ion heating rate may be similar to the electron heating rate. Thus our estimate of the ion- 
electron equipartition rate may be an underestimate by as much as an order of magnitude, 
if the wave electric field can be a factor of a few higher, and the electron Landau damping 
rate only reduced by a factor 0.5. Fuller discussion of this point is beyond the scope of this 
paper. 



4.3. Thermal Conduction 

The approach taken in this paper of following a Lagrangian plasma element out through 
the solar wind acceleration region neglects thermal conduction. Here we argue that such an 
approximation is most likely justified. The electron temperature gradient is ~ 10~ 4 K cm -1 
between 1 and 2 R Q and ~ 10~ 5 K cm -1 for heliocentric distances greater than 3R Q . The 
Spitzer-Harm heat flux is then Q ~ 10 7 ergs cm~ 2 s _1 prior to the temperature maximum (for 
an initial flow speed of 5 km s -1 ; it is lower for initial speeds of 10 and 20 km s _1 ) directed 
back towards the sun, and oppositely directed and about an order of magnitude lower once 
the wind has passed the temperature maximum. The cooling rate of the plasma occupying 
a length of order R Q between 2 and 2 R Q heliocentric distance is then ~ 0.6 s -1 , assuming 
an average electron density in this region of 10 5 cm~ 3 . This cooling rate is much faster than 
that due to adiabatic expansion, ~ 10 -4 s _1 . 

However in steep temperature gradients and collisionless plasma the Spitzer-Harm ther- 
mal conductivity is not appropriate. We estimate the collisionless heat conduction as follows: 
Begelman & Chiueh (1988) give the electron velocity diffusion coefficient in the "moderately 
non-linear regime" (i.e. corresponding to the turbulence level just before ion trapping sets 
in) as 

e 2 

D\\ «e * \ " 7 . (15) 



Putting 8E\\ = EQ\Jm e /2rrii and using equation 11 we find D\\ ~ e 2 f 2 5 2 /400m e mjC 2 c<j ~ 
4 x 10 17 (^/200kms~ 1 ) (w/4 x lO^ads" 1 ) cm s . The time taken to accelerate an electron 
from rest to a parallel velocity v\\ is t ~ v^/ D\\, which evaluates to ~ 1 s taking v\\ = 1.7 x 10 9 
cm s -1 , the electron thermal speed at T e = 10 7 K. The electron-electron collision time 
entering the Spitzer-Harm conductivity is around 3000 s, so in lower-hybrid turbulence the 
heat conduction rate should be reduced from our former estimate by a factor 1/3000, giving 
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a cooling rate 2 x 10~ 4 s _1 , which is now comparable with that due to adiabatic expansion. 
Higher degrees of lower hybrid turbulence (see above) would decrease the conductive cooling 
rate even further. 

That the heat conductivity should be reduced from the Spitzer-Harm value in steep 
temperature gradients is well known (e.g Bell, Evans & Nicholas 1981; Salem et al. 2003). 
However our estimate of heat conduction is significantly lower than one would normally 
predict. Saturated conduction is usually defined as heat conduction where electrons carry 
their own thermal energy at their thermal speed. The rationale for this is that at faster flow 
speeds (relative to the ions), a Buneman instability would develop, generating Langmuir 
waves that would inhibit the heat flow. Applied to the case above, the saturated heat flux 
would be ~ 3.5 x 10 5 ergs cm~ 2 s~ 1 . In our situation the heat conducting electrons may also 
excite lower-hybrid waves with a similar threshold drift velocity. In fact we should expect 
the limiting flux to be lower than this simple estimate, as indeed it is, because in our case 
not only heat conducting electrons but also ions generate the lower-hybrid turbulence. Any 
electron heating mechanism invoked to explain the observed charge state distributions in 
the fast solar wind must require an anomalous thermal conductivity to avoid the deposited 
heat from being conducted back to the coronal hole, which would conflict with the SUMER 
temperature diagnostics of Wilhelm et al. (1998) and David et al. (1998). We consider that 
the heating and anomalous thermal conductivity provided simultaneously by lower-hybrid 
waves to be a desirable aspect of the model, in that one mechanism provides both features. 



5. Summary and Conclusions 

In this paper we have argued that a small amount of the energy deposited in ions 
between 1.5 and 2 R Q should eventually find its way to the electrons via an instability 
that generates lower hybrid waves. No attempt has been made to address the ion cyclotron 
heating problem, other than to show that the density gradients we require are similar to 
those postulated by Markovskii (2001) to generate electrostatic ion cyclotron waves from a 
global resonant MHD mode. As such, our line of reasoning is complementary to that in a 
recent paper by Cranmer & van Ballegooijen (2003). In addressing the larger problem of ion 
cyclotron heating, these authors speculate that low frequency Alfven waves can Landau damp 
on electrons. This parallel heating should produce electron beaming and discrete phase-space 
holes which may heat ions via stochastic processes. Thus the ion heating derives from the 
electron energization (by low frequency Alfven waves), rather than the electron heating 
deriving from the ion heating as in the picture presented here. 

While a number of quantitative issues remain unresolved (atomic data, thermal conduc- 
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tion), we believe that an explanation for the observed fast wind elemental charge states in 
terms of lower hybrid wave electron heating is certainly plausible, and should be considered 
along with other possibilities already discussed in the literature (e.g Vocks & Mann 2003). 
More accurate numerical models would allow us to exploit this interpretation and allow more 
rigorous investigation of density inhomogeneities in the fast solar wind, with important con- 
sequences for the generation of ion cyclotron waves throughout the extended corona. As well 
as being potentially crucial to our understanding of the acceleration of the fast solar wind, 
the lower hybrid wave instability is also of interest elsewhere in astrophysics, and the Sun 
and solar system offer an attractive "laboratory" for the exploration of collisionless plasma 
physics processes that are otherwise inaccessible to experiment (see e.g. the recent debate 
on electron- ion equilibration in ADAFs in Binney 2003; Quataert 2003; Pariev & Blackman 
2003). 

The work was supported by NASA Contract S13783G and by the NRL/ONR Solar 
Magnetism and the Earth's Environment 6.1 Research Option. I am grateful to Steven 
Cranmer for enlightening discussions. 
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Fig. 1. — Schematic diagram illustrating the excitation of lower hybrid waves in a density 
gradient. The magnetic field direction is vertical, and density is increasing to the right. 
Gyrating ions can give a local anisotropy in the distribution function in the direction per- 
pendicular to both of these vector, into or out of the page, when the characteristic length of 
the density gradient becomes comparable to the ion gyroradius. 
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Fig. 2. — Plot of maximum lower- hybrid growth rate in units of the wave frequency 7«/u; 
(the factor Afq 2 /M is omitted) and the wavevector where this maximum is found in units 
of u/vth± against the density scale length in units of the ion gyroradius L/r g . 
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Fig. 3. — Plot of solar wind ion flow speeds adopted in the ionization models. For reference 
we also show models Bl and B2 from Cranmer et al. (1999), some flow speeds determined 
from O VI Doppler dimming observations by Patsourakos & Vial (2000), Giordano et al. 
(2000), and direct Doppler shift measurements by Hassler et al. (1999) and Wilhelm et al. 
(2000). 
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Fig. 4. — Plot of solar wind densities derived from the various models with different initial 
flow speeds, 5, 10, 20, and 60 km s -1 , increasing bottom to top. Densities derived from a 
diagnostic line ratio in Si VIII (Doschek et al. 1997) and polarization brightness (Cranmer 
et al. 1999) are given for comparison. 
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Fig. 5. — Plot of electron temperature variation with heliocentric distance for models with 
an initial flow speed of 10 km s _1 and varying degrees of collisionless ion-electron coupling, 
7' = (^jiMi/ooAfq 2 ) (cu/kviy) 2 . Measurements from the Mg IX 706/750 temperature sensitive 
ratio by Wilhelm et al. (1998) are given for comparison. We have estimated error bars on 
their points for 1.3 and 1.6i2 from the scatter of points given in their Fig. 8. 
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Fig. 6. — Plot of electron temperature variation with heliocentric distance for models with 
7' = i^/iMi/ujAfq 2 ) (u/kviy) 2 = 0.5 and initial flow speeds varying in the range 5-60 km s _1 . 
Measurements from the Mg IX 706/750 temperature sensitive ratio by Wilhelm et al. (1998) 
are given for comparison. We have estimated error bars on their points for 1.3 and 1.6i2 Q 
from the scatter of points given in their Fig. 8. 
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Fig. 7. — Plot of the evolution of the O ionization balance with heliocentric distance for 
initial flow speed 20 km s _1 and (^Mi/ 'ujAfq 2 ) (uo 2 /k 2 v 2 y ) = 0.5, corresponding to L/r g ~ 
2. The initial ionization balance corresponds to the coronal hole electron temperature of 
9 x 10 5 K. Increased ionization starts at about 1.5i? , as ion-electron energy transfer increases 
in response to the strong ion cyclotron heating at this location. Charge states are frozen in 
beyond a distance of 2-2.5 R Q , and correspond to those measured in situ by Ulysses. 



-27- 




1 2 3 4 5 6 7 

heliocentric distance (solar radii) 



Fig. 8. — Plot of the evolution of the Si ionization balance with heliocentric distance for 
initial flow speed 20 km s _1 and ('jiMi/cuAfq 2 ) (u 2 /k 2 vf y ) = 0.5, corresponding to L/r g ~ 
2. The initial ionization balance corresponds to the coronal hole electron temperature of 
9 x 10 5 K. Increased ionization starts at about 1.5i? , as ion-electron energy transfer increases 
in response to the strong ion cyclotron heating at this location. Charge states are frozen in 
beyond a distance of 2-2.5 R Q , and correspond to those measured in situ by Ulysses. 
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Fig. 9. — Plot of the evolution of the Fe ionization balance with heliocentric distance for 
initial flow speed 20 km s _1 and ('jiM i /ujAfq 2 ) [tu 2 /k 2 vf y ) = 0.5, corresponding to L/r g ~ 
2. The initial ionization balance corresponds to the coronal hole electron temperature of 
9 x 10 5 K. Increased ionization starts at about 1.5i? , as ion-electron energy transfer increases 
in response to the strong ion cyclotron heating at this location. Charge states are frozen in 
beyond a distance of 2-2.5 R Q , and correspond to those measured in situ by Ulysses. 
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Fig. 10. — Plot of maximum lower-hybrid growth rate in units of the wave frequency 
•jiMi/uAfq 2 against the density scale length in units of the ion gyroradius L/r g , for 
Maxwellian k — > oo and k — 2, 4, 6, 10, and 20. 
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Fig. 11.— Plot of V = iriiMi/uAfq 2 ) (uj/kvi y ) in units of the wave frequency against the 
density scale length in units of the ion gyroradius L/r g , for Maxwellian k — > oo and k — 2, 
4, 6, 10, and 20. 
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Table 1. C and O Ionization Fractions. 



Vstart 


7' 


c 4+ /c 5+ 


C' ,+ /C 6+ 


o 6+ /o 7+ 


km s~ 1 
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60 
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0.64 
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1.0 


0.31 


3.8 
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Geiss et al. (1995) 




0.42 


5.0 


30 


Ko et al. (1997) 
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Table 2. Mg Ionization Fractions (7' = 0.5). 



Vstart 


Mg 6 + 


Mg 7 + 


Mg 8 + 


Mg 9 + 


Mg 10 + 


km s _1 












2.5 


0.053 


0.31 


0.45 


0.14 


0.043 


5 


0.036 


0.26 


0.47 


0.18 


0.052 
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0.012 


0.16 


0.49 


0.27 


0.073 


20 


0.006 


0.11 


0.47 


0.32 


0.090 


60 


0.006 


0.10 


0.48 


0.32 


0.088 


Ko et al. (1997) 


0.028 


0.12 


0.23 


0.23 


0.40 
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Table 3. Si Ionization Fractions (7' = 0.5). 



Vstart 


Si 7 + 


Si 8 + 


Si 9 + 


Si io+ 


Si n + 


km s — 1 












2.5 


0.36 


0.43 


0.15 


0.019 


0.001 
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0.45 
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0.045 


0.003 
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0.12 


0.40 
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0.10 


0.009 


20 
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0.015 
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0.11 
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0.08 


0.31 


0.41 


0.19 


0.01 


Ko ct al. (1997) 
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0.21 
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Table 4. Fe Ionization Fractions (7' = 0.5). 



V start 


Fe 9 + 


Fe io+ 


Fe ii+ 


Fe 12+ 


Fe 13+ 


Fc 14+ 


km S _1 














2.5 


0.30 


0.34 


0.19 


0.06 


0.009 


0.001 


5 


0.14 


0.30 


0.30 


0.18 


0.049 


0.009 


10 


0.034 


0.15 


0.28 


0.31 


0.16 


0.050 


20 


0.012 


0.077 


0.21 


0.35 


0.23 


0.095 


60 


0.023 


0.11 


0.26 


0.37 


0.18 


0.051 


Geiss et al. (1995) 


0.04 


0.17 


0.28 


0.26 


0.16 


0.05 


Ko et al. (1997) 




0.16 


0.25 


0.28 


0.16 





